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ABSTRACT 

We report on the discovery of high-energy (HE; E > 0.1 GeV) and very high-energy (VHE; 
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E > 100 GeV) 7-ray emission from the high-frequency-peaked BL Lac object RBS0413. VER- 
ITAS, a ground-based 7-ray observatory, detected VHE 7 rays from RBS0413 with a statis- 
tical significance of 5.5 standard deviations (a) and a 7-ray flux of (1.5 ± 0.6 s t a t ± 0.7 sys t) x 
10~ 8 photons m~ 2 s _1 (~ 1% of the Crab Nebula flux) above 250 GeV. The observed spec- 
trum can be described by a power law with a photon index of 3.18 ± 0.68 sta t ± 0.30 sys t. Con- 
temporaneous observations with the Large Area Telescope (LAT) on the Fermi Gamma-ray 
Space Telescope detected HE 7 rays from RBS 0413 with a statistical significance of more than 
9<7, a power-law photon index of 1.57 ± 0.12 st at^o iLys an d a 7 _r ay flux between 300 MeV and 
300 GeV of (1.64 ± 0.43 st atto;22sys) x 10 ~ 5 photons uT 2 s -1 . We present the results from 
Fermi-ljAT and VERITAS, including a spectral energy distribution modeling of the 7-ray, quasi- 
simultaneous X-ray (Swift-XRT), ultraviolet (Swift-VVOT) and i?-band optical (MDM) data. 
We find that, if conditions close to equipartition are required, both the combined synchrotron 
sclf-Compton/external-Compton and the lepto-hadronic models are preferred over a pure syn- 
chrotron self-Compton model. 

Subject headings: BL Lacertae objects: individual (RBS 0413 - VER J0319+187)- 7 rays: galax- 
ies 



Introduction 



Blaz ars are active galactic n uclei that have their jet axis oriented at a small angle with respect to the 
observer ( Urrv fc Padovani 1995[) . They are observationally classified as either flat-spectrum radio quasars 
(FSRQ) or BL Lacertae (BL Lac) objects according to the broad line emission in their optical spectra. Recent 
studies interpreting the differing spectral hig h-energy (HE) 7-ray pr operties of the FSRQs and BL Lacs based 
on physical mechanisms can be found in e.g.. lGhisellini et al.l (|2009i ). Blazars are known to emit non-thermal 
radiation characterized by a double-peaked spectral energy distribution (SED). The low-energy component, 
generally covering radio to UV/X-ray bands, is usually explained as due to synchrotron emission from rela- 
tivistic electrons in the blazar jet. The origin of the HE component, occurring in the X-ray to 7-ray regime, is 
still not completely resolved and could be due to emission from a relativistic particle beam consisting of lcp- 
tons and / or hadrons. In leptonic models, very high energy (VHE) photons are produced by inverse- Compton 
(IC) scattering of low-energy photons off the synchrotron-emitting electrons. The soft seed photon s for t he 
IC process can be the synchrotron photons (synchrotron self-Compton, SS C. e.g.. iMaraschi e t al. (1992)), 



or they may originate from ambient radiation (external-Compton, EC, e.g., Dermer fc Schlickeiser 



1993)) 



Hadronic models include synchrotron emission from protons (e.g., Aharonianl J200ol)l a nd 7r°-decav from 



hadronic interactions with subsequent electromagnetic cascades (e.g.. iMiicke et aL ( 20031) ) 



RBS 0413 was discovered in the X-ray band (IE 0317.0+1834) dur ing the Einstein Medium Sensitiv- 



ity Survey and was optically identified as a BL Lac (|Gioia et al.lll984l ). The object was als o detected as 



19901 ). It exhibits sign ificant and variabl e optical polarization (|Stocke et al 



a rad io emitter with the Very Large Array of the National Rad io Astronomy Obs ervatory (jStocke et al 



Having a "feature- 



less" optical spectrum Jstocke et al. 19891) and an estimated synchrotron peak frequency logf^ Wk/Hz) 



16.99 ([Nieppola et al.l2006l ). RBS 0413 is cla ssified as a high-frequency-peaked BL Lac object (HBL. lPadovani k, Giommi 
|l995[ )). It is located at a redshift of 0.190 jGioia et al.lll984Hstocke et al.lll985[ ). 



The MAGIC Collaboration observed RBS 0413 in 2004 December-2005 February for a livetime of 6.9 

1 , at 200 GeV, assuming a power-law 



hr and reported a VHE flux upper limit of 4.2 x 10 12 erg cm 2 s 
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spectrum with a photon index of 3.0 ((Albert et al.ll2008h . VERITAS observed the source in the 2008-2009 
season and obtained a marginal significance of ~ 3a. In 2009, Jer mt-Large Area Telescope (LAT) detected 
HE emission from the direction of RBS0413 (jAbdo et al.ll2010al) . triggering new VERITAS observations. 
These new observations, combi ned with the previou s data, resulted in the detection of RBS 0413 as a VHE 
7-ray emitter in 2009 October (|Ong fc Fortinll2009h . 



2. VERITAS Observations and Analysis Results 



VERITAS is a ground-based 7-ray observatory sensitive to 7 rays with energy between 100 GeV and 
30 TeV. Located at the Fred Lawrence Whipple Observatory (FLWO) near Amado, southern Arizona, USA 
(1.3 km above sea level, N 31° 40', W 110° 57'), the array consists of four imaging atmo spheric Cherenkov 
telescopes, each having a diameter of 12m and a field of view of 3°. 5 (jHolder et alj|2008l ). During the 2009 
annual shutdown (July - August) one of the telescopes was relocated to give a more symmetrical array layout. 
In addition, a new mirror-alignment system applied in 2009 Spring contributed to an impr ovement in the 
point -spread function (PSF), with a decrease of 25% — 30% in the 80% containment radius (jMcCann et al 



2010). As a consequence, the sensitivity showed a significant improvement, and the observation time required 
for a 1% Crab Nebula detection dropped from ~ 48 hr to less than 30 hr. For observations at 70° elevation, 
the energy res olution is 15% — 20% , and the angular resolution, defined as the 68% containment radius, is 
less than 0.1° (jPerkins et al.ll2009l ). 



VERITAS observed RBS 0413 for 48 hours in total, using wobble mode (|Aharonianll2001[ ). with north, 
south, east and west wobble positions. After discarding observing runs compromised by bad weather, and a 
small number affected by hardware problems, 26 hours remained for analysis. One third of these data were 
obtained with the old array configuration (Sep 2008 - Feb 2009, MJD 54732-54883) and the rest with the 
new array (Sep 2009 - Jan 2010, MJD 55092-55485). Approximately 3 hours of data with the old array were 
taken under weak moonlight, which leads to a higher energy threshold for those observations. The source 
elevation in the data set ranges from 57° to 79°, with an average of ~ 70°. Data analysis steps consist 
of calibration, image par ameterization 



Hillas ll985T ). event reconstruction, ba ckground rejection and signal 
20081 ). For signal extraction, a 6 2 cut (jPaniel et al.ll2008h of 0.0169, 



extraction as described in lDaniel et al 
optimized for a point source of 1% strength of the Crab Nebula, was used. 



RBS 0413 is a weak source in the VHE regime. Using a "refiected-region" background estimation (jAharonian 



20011 ). an excess of 180 events and a significance of 5.5<r are obtained, for the source location at RA 



03 h 19 m 47 s ± 4| tat ± 7 s syst and deck = 18°45'.7 ± l'.0 stat ± l'.8 syst (J2000 coordinates). The VERITAS sig- 
nal is consistent with a point source, and we name the object VER J0319+187. The energy distribution 
of 7-ray events extends from ~ 250 GeV to ~ 1.0 TeV (see Table [T] for a list of spectral data points) and 
is well described by a power-law function, dN/dE = F$E~ T . The best fit is obtained with photon index 
T = 3.18±0.68 sta t ±0.30 syst and flux normalization F = (1.38 ± 0.52 sta t ± 0.60 syst ) x 10~ 7 TeV -1 xaT 2 s" 1 



at 0.3 TeV, with a value of \ 2 P er degree of freedom (x 2 /dof) of 0.14/2 (see Figure Q]). 



The integral flux above 250 GeV is (1.5 ± 0.6 8tat ± 0.7 syst ) x 10" 



corresponding to a flux level 



of approximately 1% the flux of the Crab Nebula. No significant flux variability is detected (see Figure [2] 
top panel and caption for details of t he light-curve analysis ). An upper limit (99% confidence level) on the 
fractional variability amplitude (R„ r . IVaughan et al.l (|2003l )) yields F val < 3.2. 
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Energy [TeV] 



Fig. 1. — VERITAS measured photon spectrum of RBS0413. See the text for the parameters of the power-law fit 
shown. 



Energy 


Flux 


Significance 


(TeV) 


(m-^TW- 1 ) 




0.30 


(1.3 ±0.7) x 10 -7 


2.1 


0.42 


(5.1 ± 1.9) x 10~ 8 


3.0 


0.60 


(1.5 ±0.7) x 10~ 8 


2.3 


0.85 


(4.3 ±3.4) x 10" 9 


1.2 



Table 1: Differential flux measurements of RBS0413 above 250 GeV with VERITAS. The first column shows the 
mean energies, weighted by the spectral index. The errors are statistical only. 
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3. Fermi Observations and Analysis Results 



The LAT aboard the Fermi Gamma-ray Space Telescope is a pair-convers ion 7-ray detector s ensitive 
to photons in the energy range from below 20 MeV to more than 300 GeV (jAtwood et al.l 120091 ) . The 
present analysis includes the data taken between 2008 August 4 and 2011 January 4 (MJD 54682-55565), 
which covers the entire VERITAS observation interval. Events from the Pass 6 diffuse class with energy 
between 300 MeV and 300 GeV, with zenith angle < 100°, and from a square region of side 20° centered 
on RBS 0413, were selected for this analysis. The cut at 300 MeV was used to minimize larger systematic 
errors at lower energies. The time intervals when the source was close to the Sun (MJD 54954-54974 and 
55320-55339) were excluded. The data were analyzed with the LAT Science Tools version v9r2OpC0 and 
the post-launch instrument- response functions P6_V11_DIFFUSE. The binned maximum- likelihood tools 
were us ed for significance and flux calculation ( Cash 19791 Mattox et al. 1996 ). Sources from the 1FGL 
catalog |Abdo et al.lboiObh located within a square region of side 24° centered on RBS 0413 were included 
in the model of the region. The background model includes the standard Galactic and isotropic diffuse 
emission components!!. 

A point source positi onally consistent wi th RBS 0413 is detected with a significance of more than 9a 
(test statistic, TS=89; see lMattox et al.l ([19961 )). The photon energy spectrum is best described by a power- 
law function. Replacing the power-law model with a log-parabola model does not significantly improve the 
likelihood fit. The time-averaged integral flux is 7(300 MeV < E < 300 GeV) = (1.64±0.43 sta tta22 5J J x 
10 -5 m -2 s -1 , and the spectral ind e x is 1.5 7±0.12 sta tj"n \\ ■ The spectral points were calculated using the 
procedure presented in lAbdo et al.l ( 2010a ) (see Table [2]). In the energy range 100-300 GeV, no detection 
was obtained (TS < 9) and an upper limit at the 95% confidence level was derived. Figure [2] (bottom panel) 
shows the Fermi light curve with ~ 6 month wide time bins. The upper limit point in the last time bin 
has 95% confidence level. Comparing the likelihood of a model in which the flux in each time bin is free to 
vary to one where it is assumed to be constant yields a null hypothesi s probability of 11% ; thus we find no 
evidence for variability. Details on the methodology can be found in iNolan et al.l ( 2012 ). Using the same 
method as in Section [5] and including the flux value returned by the likelihood fit in the last bin, we estimate 
-Fvar < 3.2 with 99% confidence level. 



1 http: / / fermi.gsfc.nasa.gov/ssc / data/ analysis / scitools / overview.html 
2 http: / /fermi. gsfc.nasa.gov/ssc/data/acccss/lat/BackgroundModcls. html 



Energy Flux 

(GeV) (m- 2 s" 1 TeV" 1 ) 
0.55 < 2.2 x 10" 2 

1.73 (2.6 ±0.9) x 1(T 3 

5.48 (3.9 ± 1.2) x 10 -4 

17.3 (4.6 ± 2.1) x 10" 5 

54.8 (7.7 ±4.6) x 10" 6 
173 < 4.7 x 10" 6 



Table 2: Differential flux measurements of RBS 0413 with the Fermi-L AT. The energies correspond to the bin centers. 
The errors are statistical only. 
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Fig. 2. — Top: 30-day light curve for the VERITAS data. A fit with a constant function gives a x /dof value of 
14/8, corresponding to a fit probability of 8%, consistent with the hypothesis of a constant flux. The negative flux 
point corresponding to the upper limit point in the light curve was included in the fit. Bottom: the light curve for 
the Fermi data using ~ 6 month wide time bins. The shaded areas represent the time intervals that were excluded 
to avoid solar contamination. In both graphs, the dashed lines represent the constant fit function. 
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4. Swift Observations 



The VERITAS detection triggered a Swift ([Gehrels et al.l 120041 ) target-of-opportunity observation of 
RBS 0413 on 2009 November 11, with a tot al exposure of 2 . 4 ks. All Swift-XRT data were reduced using the 
standard Swift analysis pipeline described in Burrows et al. ( 20051 ) using the HEAsoft 6.8 package. Event files 
were calibrated and cleaned following the standard filtering criteria using the xrtpipeline task and applying 
calibration files current to 2010 March. All data were taken in photon-counting mode over the energy range 
0.3-10 keV. Due to the moderate count rate of 0.3 counts s^ 1 , the data are not affected by photon pile-up 
in the core of the PSF, and partial masking of the source is not necessary. Source events were extracted 
from a circular region with a radius of 30 pixels (70". 8) centered on the source, and background events were 
extracted from a 40 pixel radius circle in a source-free region. Ancillary response files were generated using 
the xrtmkarf task, with corrections applied for the PSF losses and CCD defects. The latest response matrix 
from the XRT calibration files was applied. The extracted XRT energy spectrum was rebinned to contain a 
minimum of 20 counts in each bin. 

An absorbed power-law model, including the p/ia&Jf| model for photoelectric absorption, was fitted to 
the Swift-XRT photon spectrum. The cross-sections and abundances used the standard Xspec vl2.5 values, 
as giv en in the Xspec Analy sis ManuaQ. Using a fixed Galactic hydrogen column density, Ah = 8.91 x 10 20 
cm~ 2 (jKalberla et al.ll2005l) . the best-fit model yields a x 2 /dof value of 25.9/26. Over the energy range 0.3- 
10 keV, the best-fit photon index is T = 2.22±0.07, and the normalization at 1 keV is (33.1±2.2) keV" 1 m~ 2 
s _1 . The unabsorbed integral flux is F(0.3-10keV) = (1.69 ± 0.12) x 10~ n erg cm~ 2 s _1 in the range 0.3- 
lOkeV. The absorbed integral flux in the range 2-10keV is i^(2-10keV) = (5.81 ±0.55) x 10~ 12 erg cm~ 2 s -1 . 
No flux variability is evident over the 2.4 ks exposure. 

UVOT observations were taken in the photometric band UVM2 (2246 A) (jPoole et al.ll2008l) . The 
uvotsource tool was used to extract counts, correct for coincidence losses, apply background subtraction, 
and calculate the source flux. The standard 5" radius source a perture was used, wit h a 20" background 
region. The source fluxes were dereddened using the procedure in iRoming et all |2009l ). The measured flux 
is (2.75 ± 0.11) xl0~ 12 erg cm" 2 s _1 . 



MDM Observations 



The i?-band optical data were taken with the 1.3 m McGraw-Hill telescope at the MDM observatory on 
Kitt Peak, Arizona, between 2009 Decembe r 10 and 13 . All frames were bias corrected and flat fielded using 
standard routines in IRAE0 (jBarnesl l 19931) . and instrumental magnit udes of RBS 0413 and six comparison 
stars in the same field of view were extracted using DAOPHOT |Massev fc Davisl Il99i within IRAF. 
Physical magnitudes were computed using t he physical i?-band magnitudes of the six comparison stars from 
the NOMAD catalog (jZacharias et al.ll2005l ). assuming that the magnitudes quoted in that cat alog are exact. 



then were corrected for Galactic extinction using extinction coefficients calculated following ISchlegel et al 



3 The phabs tool applies absorption using photoelectric cross-sections. 

4 http: / /heasarc. nasa.gov/docs/software/lheasoft/xanadu/xspec/XspecManual.pdf 

5 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation. 
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(1998), taken from NEDj, and converted into v¥ v fluxes. The flux shows variations of up to ~ 30% from 
day to day, with an average of (2.47 ± 0.02) x 10~ 12 erg cm~ 2 s _1 . In the case of RBS 0413, the host galaxy 
is expected to make a substantial contribution to the observed i?-band flux. We have taken this into account 
in our SED modeling by adding a phenomenological host galaxy SED to our model. 



6. Modeling and Discussion 

The non-thermal continuum of RBS 0413 exhibits a double-peaked shape, as is typical for blazars. In this 
study, we applied three different time- independent models to the observed SED, using the contemporaneous 
X-ray, UV and optical (i?-band) data to complement the Fermi-LAT and VERITAS observations (see Fig- 
ure [3]). It should be noted that these observations were not strictly simultaneous. For all of the models, the 
emission region was assumed to be a spherical blob of size i?b, moving within the jet with a bul k Lorentz fac- 



tor r . i?b was constrained using the optical minimum variability timescale log(Ai min ) = 3.75 (jXiang fc Dai 



20071 ). where At m j n is in units of seconds. The angle between the line of sight of the observer and the jet axis, 
represented by o bs, was chosen to be equal to 1/r. This is referred to as the critical or superluminal angle, 
for which the Doppler factor equals T. The synchrotron emission was assumed to originate from relativistic 
electrons with Lorentz factors distributed between 7 m i n and 7 max , following a power law with a spectral 
index q c , under the influence of a magnetic field B. The particle-escape timescale is represented by t csc 
= r/ csc Rh/c, where rj csc is the particle-escape parameter. For each model, the parameters were adjusted to 
describe the data and achieve an equilibrium between the acceleration of the injected particles, the radiative 
cooling and the particle escape. The best-fit parameters were used to calculate the relative partition between 
the magnetic field energy density and the kinetic luminosity of relativistic particles (eBc.p = L^,/L e ^ v ) for 
each mo del. All mode l spect ra were corrected for extragalactic background light (EBL) absorption using the 



model of iFinke et al.l (|2010l ). For the optical band, a phenomenological SED reproducing the archival host 



galaxy spectral points was added to the model. 

The first model we applied assumed a pure SSC scenario. The magnetic field energy density required in 
this model is only 6% of the value corresponding to equipartition with the relativistic electron distribution 
(fBc = 0.06). The model spectrum is too hard in the Fermi band (strongly curved, with T ~ 1.5 around 
10 23 Hz) and too soft in the VERITAS band (r = 4.0), albeit within the errors in both cases. On the other 
hand, while the X-ray measurements are well reproduced, the optical (i?-band) spectrum is not. 

Next, we tested a combined SSC+EC model. The external source of photons was assumed to be an 
isotropic thermal blackbody (BB) radiation field, which may be due to a torus of warm dust with temperature 
Text = 1.5 x 10 3 K. The assumed BB infrared (IR) radiation field corresponds to a v¥ u flux of ~ 5 x 10 8 x R 2 pc 
JyHz, where R pc is the characteristic size of the IR emitter in units of parsecs. It should be noted that this 
quantity is far below the measured IR flux, thus consistent with our observations. The addition of an EC 
component improves the modeling for the optical and Fermi data compared with the pure SSC model and 
leads to values for the model parameters which are very close to equipartition (eec = 1-20). However, 
the model tends to have too sharp a cutoff in the VHE band and therefore underpredicts the VERITAS 
flux measurements. This could be remedied by choosing a much weaker magnetic field and higher electron 
energies, but the resulting system would then be very far from equipartition, with 6b c reduced by at least 
two orders of magnitude. 



http://ned.ipac.caltech.edu/ 
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Frequency (Hz) 



Fig. 3. — RBS 0413 spectral energy distribution. Absorption in the VHE region due to the EBL is taken into account 
in the fits using the model of Finke et al. ( 2010l ). The models are described in detail in the text. 
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The last model we tested is a combined lepto-hadronic jet model as described in Bottcher ( 2010l) . In this 



case, the HE component of the non-thermal emission is dominated by a combination of synchrotron radiation 
from ultrarelativistic protons (£ max > 10 19 eV) and photons from decay of neutral pions. Secondary electrons 
that are produced in various electromagnetic cascades are the origin of the low-energy synchrotron emission. 
The kinetic energy of the relativistic proton population was assumed to have a single power-law distribution 
in the energy range 1.0 x 10 3 GeV < E p < 1.6 x 10 10 GeV, with a spectral index q p = 2.4. The model is a 
good description of the overall SED, and the system is close to equipartition between the magnetic field and 
the total relativistic particle content dominated by protons (eep = 0.95). As is typical for lepto-hadronic 
models, the acceleration of protons to ultrarelativistic energies (~ 10 10 GeV) requires a high magnetic field, 
30 G in this case. Although the lepto-hadronic model provides the best description for the data, it has two 
more free parameters than the SSC+EC model and is therefore less constraining. The best-fit parameters 
adopted for all three models are summarized in Table [31 

Based on our calculations, all three models are good at describing the observed data. It appears 
that if the criterion of equipartition is taken as a reasonable measure of successful blazar emission models, 
SSC+EC is preferred over SSC fo r this HBL, which seems to be in contrast with some previous blazar 



studies. See iGhisellini et al.l ( 19981) for arguments relating the presence of an EC component with the blazar 
sequence and lAbdo et al.l |2010al) for a discussion of issues encountered in explaining blazar SEDs with a 



simple one-zone homogeneous SSC model. On the other hand, we cannot discriminate between leptonic 
and lepto-hadronic mechanisms, since the SSC+EC and lepto-hadronic models provide equally reasonable 
descriptions for the observed non-thermal continuum, and we did not detect variability in the HE and VHE 
regimes given the limited statistics. Since the synchrotron cooling timescales for electrons and protons are 
different, the detection of intraday variability would be harder to explain with a lepto-hadronic scenario and 
would accordingly favor a purely leptonic scenario. Therefore, any future observation of rapid variability 
would be helpful in distinguishing between the SSC+EC and lepto-hadronic models. 
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Parameter 


Symbol 


ssc 


SSC+EC 


Lepto-hadronic 


Electron low-energy cutoff* 




7.0 x 10 4 


5.0 x 10 4 


4.5 x 10 3 


Electron high-energy cutoff* 




1.0 x 10 6 


1.0 x 10 6 


5.0 x 10 4 


Injection electron spectral index* 


Qe 


2.4 


2.5 


2.4 


Escape-time parameter (t csc = ^csc-Rb/c)* 


f)csc 


10 


300 


100 


Magnetic field (gauss)* 


B 


0.1 


0.22 


30.0 


Bulk Lorentz factor* 


r 


20 


20 


15 


Doppler factor 


D 


20 


20 


15 


Blob radius (xlO ie cm)* 


Jib 


1.1 


1.6 


0.5 


Observing angle 


#obs 


2°. 87 


2°. 87 


3°. 82 


External radiation field E density (erg cm -3 )* 


U cx t 




6 x 10~ 7 




External radiation field BB temperature* 


Text 




1.5 x 10 3 K 




Proton spectrum low-energy cutoff (GeV)* 








1.0 x 10 3 


Proton spectrum high-energy cutoff (GeV)* 


J?p,max 






1.6 x 10 10 


Spectral index of proton distribution* 


<7p 






2.4 


Kinetic luminosity in protons (erg s -1 )* 


L P 






2.0 x 10 46 


Kinetic luminosity in electrons (erg s -1 )* 


L B (jet) 


2.97 x 10 43 


1.55 x 10 43 


6.26 x 10 40 


Magnetic field energy density (erg s _1 ) 


L B (jet) 


1.82 x 10 42 


1.86 x 10 43 


1.90 x 10 46 


Equipartition parameter 


e 


e Bc = 0.06 


£Bo = 1.20 


e Bp = 0.95 


Redshift 


2 


0.19 


0.19 


0.19 



Table 3: Best-fit parameters for the SED of RBS0413 for the SSC, SSC+EC, and lepto-hadronic models. The 
parameters that were left free are marked with an asterisk. 9 h s is the superluminal angle (see Section [6j|. 



